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PROBABILISTIC INTERPRETATION OF THE CALDERON 

PROBLEM 

PETTERI PIIROINEN AND MARTIN SIMON 


Abstract. In this paper, we use the theory of symmetric Dirichlet forms to give 
a probabilistic interpretation of Calderon’s inverse conductivity problem in terms 
of reflecting diffusion processes and their corresponding boundary trace processes. 


1. Introduction 


Electrical impedance tomography (EIT) aims to reconstruct the unknown conduc¬ 
tivity K in the conductivity equation 

(1) V • (kVm) = 0 in D 


from current and voltage measurements on the boundary of the domain D. This in¬ 
verse conductivity problem is known to be severely ill-posed, that is, its solution is 
extremely sensitive with respect to measurement and modeling errors. The unique¬ 
ness question, first proposed by Calderon has been studied in extensively. And while 
it has been answered affirmatively for isotropic conductivities in the two-dimensional 
case by Astala and Paivarinta [3], it is still unsettled, at least in full generality, in 
higher dimensions. In the previous work [22] the authors have developed a probabilis¬ 
tic interpretation of the forward problem in form of a Eeynman-Kac formula using 
reflecting diffusion processes. In this work, we extend this probabilistic interpretation 
to Calderon’s inverse problem. More precisely, we study the time-changed process X 
on dD of the reflecting diffusion process X with respect to the so-called local time on 
the boundary. For the special case of the reflecting Brownian motion on the planar 
unit disk, the Beurling-Deny decomposition of this so-called boundary process is given 
by the well-known Douglas integral 


Id 


\xn(j){x)\^ dx = 


ldDxdD\S 


{•PiO - Hv)) (47r(I - cos(^ - 77 ))) d^ dr], 


where denotes the harmonic function on D with Dirichlet boundary value (j). As 
a consequence, X is of pure jump type and its jumping mechanism, the Levy system, 
is governed by the Feller kernel ( 47 r(I — cos(^ — 77 )))“^. In this special case, X 
is the symmetric Cauchy process on the unit circle, which leads somewhat naturally 
to the following probabilistic inverse problem: Is the reflecting Brownian motion the 
unique reflecting diffusion process on the planar unit disk, whose boundary process 
is the symmetric Cauchy process on the unit circle? Clearly, this question can be 
answered affirmatively for isotropic conductivities, however, the authors are not aware 
of a probabilistic proof. 

The major goal of this work is to generalize the considerations from above in order to 
formulate probabilistic inverse problems which are equivalent to Calderon’s problem. 
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We prove that the boundary process is of pure jump type and that its jumping measure 
is governed by the Levy system of the boundary process. Moreover, we show that the 
latter is completely determined by the transition kernel density of the corresponding 
absorbing dijfusion process on D. We give explicit descriptions of both the Levy sys¬ 
tem as well as the infinitesimal generator of the boundary process. This generalizes 
results by Hsu [17] for the reflecting Brownian motion and enables the formulation of 
Calderon’s problem in the form of equivalent probabilistic inverse problems in terms 
of the excursion law of the reflecting diffusion process and the Levy system of the 
boundary process, respectively. We would like to stress the fact that in contrast to the 
reflecting Brownian motion, the reflecting diffusion process we study in this work are in 
general not in the class of solutions to stochastic differential equations. In particular, 
we can not rely on ltd calculus. 

The rest of this paper is structured as follows: We start in Section 2 by introducing 
our notation. In Section 3, we recall both the forward problem of EIT as well as its 
probabilistic interpretation from [22]. Then, in Section 4, we present two equivalent 
methods to define the boundary process of a reflecting diffusion. The first one uses 
the so-called trace Dirichlet form and its potential theory, whereas the second one 
is purely probabilistic using time change with respect to the boundary local time. 
Moreover, we prove that the Dirichlet-to-Neumann map is the infinitesimal generator 
of the boundary process. Finally, we study a Levy system of the boundary process and 
the excursions of X which leads to the formulation of three equivalent probabilistic 
inverse problems in Section 5. 


2. Notation 

Let D denote a bounded Lipschitz domain in W^, d > 2, with connected complement 
and Lipschitz parameters {tojCd), i.e., for every x € dD we have after rotation and 
translation that dD n B(x^ro) is the graph of a Lipschitz function in the first d — 1 
coordinates with Lipschitz constant no larger than cd and D (1 B(x,ro) lies above the 
graph of this function. Moreover, we set Ml := {a: G M'’* : a; • < 0}, with v = Cd the 

outward unit normal on where we identify the boundary of Ml with with 

straightforward abuse of notation. 

For Lipschitz domains, there exists a unique outward unit normal vector v a.e. on 
dD so that the real Lebesgue spaces LP{D) and LP{dD) can be defined in the standard 
manner with the usual norms IJ-Up, p = 1,2, oo. The standard inner-products 
are denoted by (•, •) and (•, ■)gD, respectively. The d-dimensional Lebesgue measure is 
denoted by m, the (d — l)-dimensional Lebesgue surface measure is denoted by a and 
j • j denotes the Euclidean norm on M"^. 

All functions in this work will be real-valued and derivatives are understood in 
distributional sense. We use a diamond subscript to denote subspaces of the standard 
Sobolev spaces containing functions with vanishing mean and interpret integrals over 
dD as dual evaluations with a constant function, if necessary. For example, we will 
frequently use the spaces 

Ht^/\dD) := G = o} 

and 

HliD) :={f,€H\D):{cf,l)=0}. 

Moreover, we will frequently assume that dD is partitioned into two disjoint parts, 
diD and d 2 D. We denote by Hq{D U diD) the closure of Cf°{D U diD), the linear 
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subspace of C°°{D) consisting of functions (p such that supp((^) is a compact subset 
oi D U diD, in H^{D). For the reason of notational compactness, we use the Iverson 
brackets: Let 5" be a mathematical statement, then 


[5] 


1, if /S' is true 

0 , otherwise. 


We also use the Iverson brackets [x € B] to denote the indicator function of a set B, 
which we abbreviate by [B] if there is no danger of confusion. 

In what follows, all unimportant constants are denoted c, sometimes with additional 
subscripts, and they may vary from line to line. 


3. The forward problem and its probabilistic interpretation 


We assume that the, possibly anisotropic, conductivity is defined by a symmetric, 
matrix-valued function k : D ^ with components in L°°{D) such that k is 

uniformly bounded and uniformly elliptic, i.e., there exists some constant c > 0 such 
that 


(2) c ^ ■ k{x)^ < c|.^p, for every ^ € R"* and a.e. x £ D. 

The forward problem of electrical impedance tomography can be described by differ¬ 
ent measurement models. In the so-called continuum models the conductivity equation 
( 1 ) is equipped with a co-normal boundary condition 

(3) := KP • Vu|aD = / on dD, 


where / is a measurable function modeling the signed density of the outgoing current. 
The boundary value problem (1), (3) has a solution if and only if 

(4) (/; l)a_D = 0- 

Physically speaking, this means that the current must be conserved. Given an appro¬ 
priate function /, the solution to (1), (3) is unique up to an additive constant, which 

_ ^ /2 

physically corresponds to the choice of the ground level of the potential. If / G , 

then there exists a unique equivalence class of functions u £ 7L^(£))/R that satisfies 
the weak formulation of the boundary value problem 


/ K\Ju-\Jvdx = {f,v\QD)QD 
Jd 


for all V £ H^{D)/R, 


where v\dD '■= "fv and 7 : H^{D)/M. —>■ /M. = (iJo is the standard 

trace operator. Note that we occasionally write v instead of v\dD for the sake of 
readability. 

In his seminal paper [ 8 ], Fukushima established a one-to-one correspondence be¬ 
tween regular symmetric Dirichlet forms and symmetric Hunt processes, which is the 
foundation for the construction of stochastic processes via Dirichlet form techniques. 
Therefore we assume that the reader is familiar with the theory of symmetric Dirichlet 
forms, as elaborated for instance in the monograph [ 10 ]. 

Let us consider the following symmetric bilinear forms on L^(D): 



kS/v(x) ■ 'S/w{x) dec. 


v,w£ V{£) := H^{D) 


(5) 


£{v,w) : 
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and for the particular case k = which is of special importance, we set 

(6) £^^{v,w) ■.= - f Vv{x) ■'Vw{x) dx, v,w € V{£^^) := H^{D). 

2 Jd 

The pair {£, £>{£)) defined by (5) is a strongly local, regular symmetric Dirichlet form 
on L‘^{D). In particular, there exist an f-exceptional set Af C D and a conservative 
diffusion process X = > 0},Pa;), starting from x G D\J\f such that X is 

associated with {£,!){£)). Without loss of generality let us assume that X is defined 
on the canonical sample space = C'([0, oo); D). It is well-known that the symmetric 
Hunt process associated with ( 6 ) is the reflecting Brownian motion. Therefore, we call 
the symmetric Hunt process associated with (5) a reflecting diffusion process. 

Let us briefly recall the concept of the boundary local time of reflecting diffusion 
processes, see, e.g., [4, 15, 22]. If the diffusion process is the solution to a stochastic 
differential equation, say the reflecting Brownian motion, then the boundary local 
time is given by the one-dimensional process L in the Skorohod decomposition, which 
prevents the sample paths from leaving D, i.e., 

1 

(7) Xt = x + Wt- - iyiXs)dLs, 

2 Jo 

P 2 ,-a.s. for q.e. x G D. This boundary local time is a continuous non-decreasing process 
which increases only when Xt G dD, namely for alH > 0 and q.e. x G D 

Lt= [ [dD]{Xs)dLs, 

Jo 

Pa,-a.s. and 

E, [ [az?](X,)ds = 0 . 

Jo 

Although the reflecting diffusion process associated with (5) does in general not ad¬ 
mit a Skorohod decomposition of the form (7), we may still define a continuous one¬ 
dimensional process with these properties. More precisely, by the Lipschitz property 
of dD, we have that Dr\B{x,rD) = {{x,Xd) '■ Xd > 'y{x)}r\B{x,rD) and the Lipschitz 
function 7 is differentiable a.e. with a bounded gradient. In particular, we have for 
every Borel set B C dD fl i?(a;, r^) that 

(j{B)= f (1-f |V7(i)p') ^ dx 

d{£:($,7(£))eB} ' 

and a straightforward computation yields that the Lebesgue surface measure cr is a 
smooth measure with respect to {£,D(£)) having hnite energy, i.e.. 


juj d(T(a;) < c||u||£j for all v G 'D{£) fl C{D), 


IdD 


where we have used the inner product fi(-, •) := £{■, ■) -I- (•, •). 


Definition 3.1. The positive continuous additive functional of X whose Revuz mea¬ 
sure is given by the Lebesgue surface measure cr on dD, i.e., the unique L G such 
that 

( 8 ) lim y [ Ea;| [ (l){Xs)dLs\4’ix)dx = [ (l)ix)'ip{x) da{x) 

i-s-0+ t Jl Jq ) Jqjj 

for all non-negative Borel functions (j) and all a-excessive functions ip, is called the 
boundary local time of the reflecting diffusion process X. 
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It has been shown in the recent work [22] that the ^-exceptional set J\f is actually 
empty. 

Proposition 1 ([22, Proposition 1]). p € C°’'^((0,T] x D x D) for some 5 € (0,1), 
i.e., for each fixed 0 < to <T, there exists a positive constant c such that 

(9) \p{,t 2 ,X 2 ,y 2 ) - p{ti,xi,yi)\ < c{y/t 2 - h + \x 2 - sij + \y 2 - yi\Y 

for all to < ti < t 2 < T and all {xi,yi)^ {x 2 ,y 2 ) € D x D. Moreover, the mapping 
1 1 —>■ p{t, •, •) is analytic from (0, oo) to x D). 

By [9, Theorem 2], the existence of a Holder continuous transition kernel density 
ensures that we may refine the process X to start from every x € D hy identifying the 
strongly continuous semigroup {Tt,t > 0} with the transition semigroup {Pt,t > 0}. 
In particular, if v is continuous and locally in the Fukushima decomposition 

holds for every x € D, i.e., 

(10) v{Xt) = v(Xo) + Mf + Nf, for all t > 0, 

Pa;-a.s., where M"" is a martingale additive functional of X having finite energy and 
N'’ is a continuous additive functional of X having zero energy. 

Moreover, both M'’ and N'" can be taken to be additive functionals of X in the 
strict sense, cf. [10, Theorem 5.2.5]. 

Finally, note that the 1-potential of the Lebesgue surface measure a of dD is the 
solution to an elliptic boundary value problem on a Lipschitz domain with bounded 
data. By elliptic regularity theory, cf., e.g., [11], this solution is continuous, implying 
that the boundary local time L exists as a positive continuous additive functional in 
the strict sense, cf. [10, Theorem 5.1.6]. 

For the probabilistic interpretation of the Neumann boundary value problem corre¬ 
sponding to the continuum model, we require the following assumption on the conduc¬ 
tivity k: 

(Al) There exists a neighborhood U of the boundary dD such that k\i( is isotropic 
and equal to 1. 

(A2) There exists a neighborhood U of the boundary dD such that k\u is Holder 
continuous. 

Remark 1. Notice that, theoretically, the assumption (Al) imposes no restriction to 
generality. More precisely, it can be shown using extension techniques that for domains 
D,D C such that D C D, the knowledge of both, the Dirichlet-to-Neumann map 
Afi on dD and yields the Dirichlet-to-Neumann map A^ on dD. 

The main result for the continuum model (1), (3) is the following theorem from [22]. 

Theorem 3.2 ( [22, Theorem 6.4.]). Assume that n satisfies (Al) or (A2) and let 
D denote a bounded Lipschitz domain. Let f be a bounded Borel function satisfying 
(/, 1)qjj = 0. Then there is a unique weak solution u € C{D) r\Hl{D) to the boundary 
value problem (1), (3). This solution admits the Feynman-Kac representation 

(11) u{x) = lim Ea; [ f{Xs)dLs for all x € D. 

Remark 2. We would like to emphasize that the regularization technique employed 
in the proof of [22, Theorem 6.4.] may be easily modified to prove the Feynman-Kac 
formula 


u{x) = Ex(l){Xr(D)), X e D 
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for the conductivity equation (1) with Dirichlet boundary condition u\dD = 4>i where 
(j) G and 

t{D) := inf{t > 0 : e R‘^\D} 

denotes the first exit time from the domain D. This follows from the Lipschitz property 
of dD, implying that all points of dD are regular in the sense of [18, Chapter 4.2]. 


4. Boundary processes of reflecting diffusions 


4.1. Definition and properties. As we are going to use regularity results that are 
not readily available for general Lipschitz domains, we assume throughout this chapter 
that D has a smooth boundary in order to avoid technical difficulties. However, we 
expect our results to hold for general Lipschitz domains. 

Let denote the absorbing diffusion process on D which is obtained from the 
reflecting diffusion process A on D by killing upon hitting of i.e.. 


( 12 ) 




At, t< t{D) 
d, t > t{D). 


By the Markov property of A, possesses a transition kernel density which may be 
expressed as 

(13) p°(t,x,y) =p(t,x,y) - E^{p(t - t{D), Xr(D),y)[TiD) < t]} 

and the regular symmetric Dirichlet form on Lf{D) associated with A° is (f,iLg(Zl)), 
cf. [10]. 


Lemma 4.1. Let k : D ^ be a symmetric, uniformly bounded and uniformly 

elliptic conductivity with components Kij € i,j = 1, ...,d, such that k satisfies 

(Al). Then the absorbing dijfusion process A® possesses a transition kernel density p'^ 
with the following properties: 

{i) is jointly Holder continuous with respect to (t,x,y); 

(ii) p^ is in Hq{D) as a function of x and y, respectively; 

(iii) p'^\u € C^(IL U dD) fl C^fU) as a function of x and y, respectively, where U 
denotes the neighborhood of dD from assumption (Al); 

(iv) p^ is continuously differentiable with respect to t as a Banach space valued map. 


Proof. The property (i) follows directly from Proposition 1 by (13). 

To show the property {ii), first note that for every x G D, we have Vx,t{y) = 
p^{t/2,y,x) is Holder continuous for every t > 0 and in particular it is bounded and 
in Lf{D) since D is bounded. This implies that for every t > 0 we know that 

Tt/ 2 Vx{y) = SyVx{Xt/2)[t/2 <td]= [ p°{t/2,y,z)p°{t/2,y,x)dz = Vx, 2 t{y) 

JD 

and moreover, Tt/ 2 Vx G D{C) = H({D). This implies the claim. 

The property (in) is a direct consequence of assumption (Al) and the regularity of 
the fundamental solution for the heat equation, cf. [15]. 

Finally, (iv) follows by the same reasoning we used in the proof of Proposition 1, 
see [22]. □ 


A function u G Hi(^(D) is said to be C-harmonic, provided 
(14) Ea;|'u(A^^gj)j < oo and u{x) = for all x G D 
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and every relatively compact open subset D C D, which is clearly equivalent to 

(15) I KVu-Vuda; = 0 for all u € C'c(-D), 

Jd 

cf. Remark 2. By the same remark, the C-harmonic extension operator 'H : L°°{dD)r\ 
^i/ 2 -£(^£)) eg [0,1/2), dehned by 

n4>{x) xeD 

is well-defined. If e g (0,1/2), this is due to the fact that the Dirichlet problem 

(16) V • (kVu) =0 in H, u\dD = 

admits a unique solution u g satisfying (15), which follows from a 

saddle point formulation introduced in [20]. 

The following lemma generalizes a result from Aizemnan and Simon [1] for the 
absorbing Brownian motion. 


Lemma 4.2. Let n : D ^ be a symmetric, uniformly bounded and uniformly 

elliptic conductivity and suppose, in addition, k satisfies (A2). Then for every bounded 
Borel function (j) on dD we have 

E,,(l){Xr(^D))[T{D) <t]= [ f (j){y)d {y)P°{s,x,y)<lsda{y), x € D. 

JdD Jo 

Proof. We will proceed in the spirit of the Aizenman and Simon [1] proof but we will 
reformulate it to cover our case. Let g (^/“([O,^)) be a smooth approximation of 
the indicator function [0 < s < t] and vo{y) = p^{e,x,y) for given x € D and e > 0. 
Note that vq G L‘^{D). 

Let w g H^{D). Then since Vs = TgVo satisfies an abstract Cauchy problem, cf., 
e.g., [21], we have by the equivalent variational formulation that 

poo poo 

/ {Cvs,w)'tlj{s)ds = / (Us,w)'!/(s)ds-|-Ea,w(Ae)[e < t] 

Jo Jo 

Let ipe be a non-negative Hq (D) function such that is an approximation of the indicator 
function [Hg] or more precisely ipe[D 2 e] = [D 2 e] and (pe[D \ De] = 0. Rewriting w on 
the right-hand side by tc = + wi^ — Te) we get by the definition of weak derivative 

that 

{Cvs,w) = -(kVus, V((/?eU;)) -f {Cvs, (1 - (pe)w) 

= -{kXVs,Xw) -I- (kVUs, V((l - <Pe)w)) + {CVs, (1 - (Pe)w) 


If w is an £-harmonic, the first term on the right is 0. If k is Holder continuous close 
to the boundary dD, then the elliptic regularity shows that 

lim|(£us, (I - ((£)e)w)l = 0 
€—>•0 

since (I — ipe)w —>■ 0 in H^~^{D). The term 

(kVUs, V((I - ipe)'w)} = (kVVs, (1 - ipe)Vw)} - (kVVs , wV(pe} ■ 

Since w g H^(D), we have 

limK^Vus, (1 — (/Je)VM;)] = 0 

since (1 — ipfjXw —>■ 0 in Lf{D). Since the boundary dD is Lipschitz and compact, it 
can be further divided into a finite partition of H^{D) and within this partition the 


Petteri Piiroinen and Martin Simon 


V(^e{x) = —e ^[x £ De \ D 2 e]i'{x*) where x* £ dD is the unique point such that after 
the Lipschitz change of coordinates, x = ^^{x*). Therefore, 

- lim|(KVus, wV(^e)| = / d^^ty)Vs{y)w{y)da{y) 

<^^0 JdD 

and all in all, we have shown that 

pOO p pOO 

/ tp{s) d^^i^y)Vs{y)w{y)dcF{y)ds = I (u*,■u;)i/i(s) ds + Ea;w(X<;)[e < r ] 

Jo JdD Jo 

+ 0 ( 1 ). 

Since the semigroup {T^, t > 0} is continuous and ip is an approximation of the indicator 
function [0 < s < t], we further deduce that 

rt 


n dK^{y)Vs{y)w{y) do-(y) ds = -(ut, w) + w{x) + o(l) 

>D 

D. Since 

{vt,w)= vt{y)w{y)dy= / p°{t,y, z)p°{e,x, z)w{y) dy dz = Tt+ew{x) 
J D Jd Jd 


Jo JdD 
for m-a.e. x £ D. Since 


Jd 

we have deduced that 

/ / di^^(y)Vs{y)w{y) da{y) ds = -Ttw{x) + w{x) + o(l) 

Jo JdD 

for m-a.e. x £ D and for every £-harmonic function w £ H^{D). Using the Chapman- 
Kolmogorov and the analyticity of p^ with respect to t we can write this as 

/ / dKv{y)P°{s, X, y)w{y) dcr{y) ds = -Ttw{x) + w{x) + o(l) 

J e J dD 

which implies that 

(17) / / d^^(^y-)P°{s, X, y)w{y) da{y) ds = -Ttw(x) + w{x) 

Jo JdD 

for m-a.e. x £ D. We can use the Markov property (as in the proof of of [22, 
Lemma 6.3]) to show that identity (17) holds for every x £ D and for every £-harmonic 
function w £ H^{D). In particular, we can reformulate this as 

/ / d^^(^y)p°{s, X, y)4i(ij) duijj) ds = -TtUcpix) + Hcpix) 

Jo JdD 

= -E:cn(j){Xt)[t < t] d-Hcpix) 

for every (p £ H^^‘^{dD). By Markov property and the fact that for every x £ D the 
stopping time r > 0 a.s. we have 

E^'HcP{Xt))[t <t]= E^E{cP{Xr)[t < r] I J-*) = E,<^(X,)[t < rj. 

Thus, the right-hand side of the identity (18) is by the definition of the harmonic 
extension operator 

-Ej;(p{Xr)[t < t] -\-E^(p{Xr) = E^(p{Xr)[t > t]. 

We have therefore shown that 

(19) f f d^^^y)p°{s,x,y)(p{y)da{y)ds = Ea;(p{Xr)[T <t] 

Jo JdD 


(18) 
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for every cj) € H^/^{dD). This implies by the Monotone Class Theorem that the 
identity (19) holds for every bounded and measurable function (j). □ 

Lemma 4.2 yields the joint distribution of the pair with respect to 

the measure P^;, a; £ D, namely 

(20) f‘^{T{D) £ At, £ dy} = x, y) At Aa{y). 

In particular, there exists a Poisson kernel with respect to the Lebesgue surface mea¬ 
sure (T which is given by 

nOO 

K^{x,y):= / x € D, 

Jo 

such that the hitting probability of i? £ B{dD) for X starting in a; £ D is 

^x{Xr(D) & B} = [ K^{x,y)da{y). 

JB 

Following [6], the so-called trace Dirichlet form {£, !){£)) of {£,V{£)) is given by 

(21) £{v,w) :=£{nv,nw), V{£) ■=Ve{£)\aD n L^{dD). 

It is a symmetric regular Dirichlet form on Lf{dD), cf. [10, Theorem 6.2.1], whose 
domain is characterized by the following lemma. 

Lemma 4.3. Let k : D ^ be a symmetric, uniformly bounded and uniformly 

elliptic conductivity. Then'D{£) = H^^^(dD). 

Proof. By the standard trace theorem it suffices to show 'De{£) C P[^{D) = 'D{£). 
Let V £ 'De{£)- Therefore, by the definition of the extended Dirichlet space, there is 
a sequence {vk}ken C H^{D) such that (vk) converges m-a.e. on D to u and (vk) is 
f-Cauchy sequence. 

Since (vk) is a f-Cauchy, we have by the Poincare inequality and the boundedness 
of K that 

(22) II Un Vjn AlVn A (D) — ^11^ '^m)||2 — ^l£i'^n Xjri), 

where Mu := \D\~^{u, 1)2[L*]- In other words, {z;„ —is an iJ^(D)-Cauchy se¬ 
quence and thus exists a tc £ H^{D) such that Vn — XIvn —?► rc in H^{D). Furthermore, 
we may choose a subsequence {wn — Mwn} C {u„ — Mvn} such that Wn — Mwn —>■ w 
for m-a.e. in D. On the other hand, we have that Wn —^ v for m-a.e. in D and thus, 
Xlwn = Wn — (wn — Mwn) w — V ior m-a.e. in D. Since Xlwn is a constant function 
for every n, we deduce that w — v = c[D]. Therefore, v = w + c[D] £ H^{D) as 
claimed. □ 

Now let us consider the Hunt process X = (D, T, {X, t > 0}, Px) associated with the 
trace Dirichlet form, which is a cr-symmetric pure jump process on dD. More precisely, 
as the energy measure of M'', v £ 'De{£) = H^{D) does not charge the boundary dD 
and due to the existence of the Poisson kernel, we may apply [6, Theorem 6.2] to obtain 
for v,w £ V{£) n Cc{dD) the Beurling-Deny decomposition 

(23) £(v,w)=( {v{x)-v{y)){w{x)-w{y))AJ{x,y). 

JdDxdD\5 

The jumping measure dff (x, y) is determined by the absorbing diffusion process Xq 
and can be characterized by the so-called Feller measure which depends only on the 
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associated symmetric strongly continuous contraction semigroup on Lp‘{D) and the 
Poisson kernel, cf. [6]. X will be called the boundary process of X. 

Note that the boundary process of X may be equivalently constructed using the gen¬ 
eral theory of time changes of diffusion processes with respect to positive continuous 
additive functionals: We have seen that the boundary local time L is a nondecreasing, 
{Xut > 0}-adapted process that increases only when X is at the boundary. Follow¬ 
ing [10], we define the right-continuous right-inverse r of L by 

(24) r(s) := sup{r > 0 : < s}- 

The random variable t ( s ), s € [0, oo), is a stopping time with respect to the right- 
continuous history > 0} of X since { t ( s ) > t} = {Lt < s} G Tt and moreover, 

by continuity of X we see that for every s G [0,oo) the process X is at the boundary 
dD at time t(s). Therefore, we can equivalently define the boundary process X of X 
as the time-changed trace 

:= X,(t) 

and the boundary filtration 

Tt := Trit)- 

To be precise, we know that the boundary local time L is a positive continuous 
additive functional in the strict sense since the corresponding 1-potential is bounded 
by elliptic regularity and therefore, the boundary dD coincides with the so-called quasi- 
support of L, cf. [10, Theorem 5.1.5]. Moreover, since dD is smooth, every boundary 
point is a regular point, cf. [18], so that the boundary process is a cr-symmetric Hunt 
process on the boundary dD^ cf. [10, Theorem A.2.12., Theorem 6.2.1]. 

Remark 3. Note that due to the refinement obtained obtained by Proposition 1 for 
the reflecting diffusion process X, the boundary process X can be defined without 
exceptional set, i.e., for every starting point x G dD. 

We note that the representation Theorem 3.2 can be expressed with the help of the 
boundary process X, the first exit time td and the first exit place Xt^. 

Lemma 4.4. We have that 

pT{b) pb 

(25) / f{s)dLs= /(T(s))ds 

JT(a) Ja 

for every bounded and measurable function f: R —^ R. 

Proof. This follows by Monotone Class Theorem from the observation that 

[Lai L^,] — TO Pa^b: 

where we have set ga,b{t) ■= [t G [a,5]]. □ 

4.2. Infinitesimal generator. The following theorem yields the infinitesimal gener¬ 
ator of the boundary process. 

Theorem 4.5. Let k : D ^ be a symmetric, uniformly bounded and uniformly el¬ 

liptic conductivity. Then the infinitesimal generator of X is the Dirichlet-to-Neumann 
map A^. 
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Proof. We will give two different proofs for this result. One is with the trace Dirichlet 
forms and one with the help of Theorem 3.2 and Change of Variables lemma 4.4. 

Let v,w G n(£i) = H^/'^{dD). Therefore, 

£{v,w) = £{'Hv,'Hw) = I w{x){dKi,(^j;'j'Hv{x))\dD d<7{x) 


IdD 


w{x){d^^,(^^'j'Hv{x))\gD d(j{x) 
an 

= {w, h.K,v) 

We can factorize A„ = A*A, where A: H^/‘^{dD) —>■ L^{dD). The A can be seen as an 
unbounded operator in Lf{dD) with domain £>{£)■ Thus, 

£{v,w) = {Aw,Av)L'^(dD), 

which implies A = 'J—C or = —£, where C is the infinitesimal generator of X 
(c.f. [10]). 

The another proof is even more probabilistic in nature. Let G C°°{dD) and 
V = AkV^. We know that v is Holder continuous and {v, l)dD = 0. Therefore, we may 
define 


9r^^') '•— ^xAjl — Ea; 


u(Vs)ds and goo{x) = lim gR{x). 

R^oo 


Note that the Theorem 3.2 together with the Change of Variables Lemma 4.4 gives 
goo = A;7^u = ip and in particular, g^o £ C°°{dD) C '£>{£). The transition operator Tt 
applied to g gives by Markov property 

TtgR{x) = = Ea(AK+t — At). 

Therefore, 

(Tt — I)gR{x) = —EaAt + Ea(Afl-|_t — Ar) 
and the Theorem 3.2 together with dominated convergence theorem imply that 

{Tt - I)p{x) = -EaAt = - [ Tsv{x) ds. 


This in turn implies 

C(f{x) = -v{x) = -Ak.p{x) 
which implies the claim since {£,£>{£)) is regular. 


□ 


As in [17], we show that the Dirichlet-to-Neumann map is an integro-differential 
operator in the sense of Lepeltier and Marchal [19] . 

Theorem 4.6. Let k : D ^ be a symmetric, uniformly bounded and uniformly 

elliptic conductivity with components Rij £ C^'^{D), i,j = 1, ...,d such that k satisfies 
(Al). Then the Dirichlet-to-Neumann map A^ is of form 

= A„id • Vt(/> + Af^cj) for all (j) £ H^(^{dD), 
where is given by the integral operator 

A,^cj){x) := I {cj}{y) - cj}{x) - Vrcfix) ■ {y - x)) y) d<7{y) 

JdD 

for a.e. x € dD and 


Nn{x,y) := d,^^t^,,)K^,{x,y), x,y € dD. 
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Proof. By a density argument, we may assume that G C°°{dD) fl H^^'^{dD). The 
unique solution in H^{D) of the Dirichlet problem (16) is by Lemma 4.2 

(f{y)9Ku{y)P°{t, X, y) dt dcr(y). 

Therefore, maps (j) to 

(26) A^(j){x) = 

Let us extend (j) into the neighborhood of the boundary as constant along the conormal 
direction. Therefore, the first order tangential derivative V := VT<t> gets extended 
simultaneously. We will denote the extensions (j) and V, respectively. We compute the 
conormal derivative of the function 

d 

w — ucf-mi- 

i=i 

where {W} is a vector field on the boundary defined by W{y) := yr as the projection to 
the tangent plane going through the point y. By construction, the conormal derivative 
commutes with multiplication by the extended functions and vector fields. Therefore, 

d 

= K^4> - = K,4>-V- A„W, 

1=1 

where A„1 = 0 since u{x) = 1 in H is the unique solution to the corresponding Dirichlet 
problem and therefore, the conormal derivative vanishes identically. As we have shown 
the existence of a Poisson kernel, we can write the left-hand side in a different way, 
namely 

S/w{x) = S/x f {cl){y)-^{x)-V{x) ■ {W{y)-W{x)))Kf,{x,y)da{y) 

JdD 

for almost every a; in a neighborhood of the boundary. 

Next, we show that the function y i-i Nk.{x, y){\y — x\^ A1) is integrable with respect 
to the surface measure a for every x £ dD. When k = 1, the proof given in [17] yields 
the claim for the kernel A^i corresponding to /t = 1. As we assume that k satisfies (Al), 
the operator A^ — Ai is a smoothing operator which follows by the standard elliptic 
regularity, cf. [14]. This implies that the kernels Ni and N have the same leading 
singularities and the claim thus follows from the estimate for Ni. As a consequence, 
we can use the dominated convergence theorem to take the differentiation inside the 
integration and we obtain 

dK.vw{x)= / {4>[y) - (i){x)-V t4>{.x) ■ {y - x))Ni.,{x,y)da{y) = Ai.,(l){x) 

JdD 

for a.e. x £ dD. □ 

5. Novel probabilistic formulations of the Calderon problem 

By the considerations from above, the boundary process is uniquely determined by 
the absorbing diffusion process so that Theorem 4.5 leads to the following proba¬ 
bilistic interpretation of Calderon’s problem: Given a boundary process X associated 
with the regular asymmetric Dirichlet form {£,H^/‘^{dD)), is X'^ the unique absorb¬ 
ing diffusion process on D such that X is the boundary process of the corresponding 
reflecting diffusion process X on D? 


^.(^{x) =¥.x4>{Xrr,) = [ [ 

JdD Jo 


Probabilistic Interpretation of the Calderon Problem 


13 


The Calderon problem in 2 dimensions is known to be solvable for isotropic k S 
L°°{D). Given the boundary process X = we can thus uniquely determine the 
generator A = A^o- The celebrated result of Astala and Paivarinta [3] says that 
whenever A„ = A and both k and kq are isotropic, uniformly bounded and uniformly 
elliptic, then k = kq. Therefore, the equality X^ = X^g must hold as well. 

The recent result by Haberman and Tataru [13] implies the same for three and 
higher when k and kq are assumed to be C^{D) or if they are Lipschitz continuous 
and close to identity in certain sense. In recent preprints, Haberman [12] improved 
this to the even conductivities for the dimensions three and four, while Caro and 
Rogers [5] improved the Haberman and Tataru technique to prove the uniqueness for 
the Lipschitz case in general in any dimension. 

When the conductivity is not assumed to be isotropic the uniqueness has always 
an obstruction, namely we have A^o = A^j, whenever ki = is the push-forward 
conductivity by a diffeomorphism F on D that leaves the boundary dD invariant. In 
the plane, this is known to be the only obstruction by the result of Astala, Lassas and 
Paivarinta [2] which holds without additional regularity assumptions on the conduc¬ 
tivity. In higher dimensions the question is still very much open in general, see [7] for 
further discussion. 

It should be emphasized that these results from analysis all rely on so-called complex 
geometric optics solutions and the authors are not aware of any probabilistic interpre¬ 
tation of such solutions. Therefore, it is an open problem to find a probabilistic solution 
to the probabilistic formulation of Calderon’s problem. 

Let us elaborate a bit on this problem by providing three seemingly different but 
equivalent versions. 

5.1. First version. It is an immediate consequence of Theorem 4.6 that the jumping 
measure of X can be described with the help of [19, Theoreme 10]. 

Lemma 5.1. Let k : D ^ be a symmetric, uniformly bounded and uniformly 

elliptic conductivity with components Kij € i,j = 1, ...,d, such that k satisfies 

(Al). Then for every non-negative Borel function cj) : dD x dD —>■ R+, vanishing on 
the diagonal, and any stopping time f of X, we have 

(27) E, V(/.(X,_,X,)[X,_ ^X«] =E, / [ (f(Xs,y)N^iXs,y)daiy)ds. 

s<f 

Proof. Suppose first that diam(II) < 1. We note that the operator in Theorem 4.6 
coincides with the integral operator causing the jumps, li tp G C°°{dD) fl H^^'^{dD) 
and it is continued as ^ S so that ip and its tangential derivative are continued 

as constants along the conormal directions in the neighborhood of the boundary dD, 
then for every x G dD, we have 

= [ {^^{x + z) - p>{x) - [\z\ < l]z ■y'iij{x))SK{x,dz), 

where we have set SK,{x,dz) := Nk,{x,x z)dax{z) and ax{B) := a{x B) for every 
Borel set B G In the same way we can extend the drift term so that the cor¬ 
responding integro-differential operator is the infinitesimal generator of an extended 
process X on the whole space K"’*. Since we know that Xt G dD for all t > 0, it follows 
that the extended process X will stay on the boundary dD if we start it from the 
boundary and it coincides with X there. 
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For this extended process X we can apply the result [19, Theoreme 10] and we 
obtain 

E,J^^(Xs-,Xs)[Xs- 

S<T 

= E^ [ [ (j){Xs,Xs + y)Ni^{Xs,Xs + y)d(Tj^ {y)ds 

Jo dR<l\{0} 

for any non-negative Borel function cj) : dD x dD —> IR+ vanishing on the diagonal and 
any stopping time f oi X. The claim follows now in this special case by change of 
integration variable. 

The general case follows by scaling: Let us denote X^ := R~^Xt. This is a reflecting 
diffusion process corresponding to on a domain D^, where := R~^D and 
K^{x) := R~^k{Rx). Since the diameter of is one, the claim holds for the boundary 
process X^ of X^. 

Let denote the local time of X^ on the boundary dD^. By definition, this is in 
Revuz correspondence with the surface measure of the boundary dD^. By using 
the Revuz correspondence and change of variables, it follows that 

Lf = RLf 

Therefore, the right-inverse of the local time has a scaling law 

r«(f) = 

which in turn implies that the boundary processes scale by the law 

vR D—1 V 

Xt — H 

and that y is an X-stopping time if and only if Ry is an X^-stopping time. 

If we denote by the conormal derivative of the Poisson kernel of X^ multiplied 
by 2 and compute the scaling law, we find out that 

= R‘^~^N^{Rx,Ry). 

With all these scaling laws, we are now ready to prove the claim for X. Let y) := 

(j){x, y) and let y be an X-stopping time. We have 

E,^0(X„X,_)[X, ^ X,_] = ^ 0«(Xf,Xf_)[Xf ^ Xf_], 

s<77 s<r]R 

where denotes the expectation given X^ = R~^x. By the first part of the proof, 

the right-hand side is equal to 

En-^J / ^^iX^,y)N^iX^,y)daiy)ds. 

Jo JdD^ 

With the change of variables y' = Ry and s' = sR~^ and the scaling law R~^y) 

R‘^~^Ni^{Xs,y), the claim follows. □ 

This result states that the pair {Nk,{x, y) dcr(t/), idt) is a Levy system of the Hunt pro¬ 
cess X. Since the positive continuous additive functional idt = t with respect to X has 
the Revuz measure a, we obtain from [6, Theorem A.l (iii)] that y) dcr(y) d(7(a:) 

coincides with the jumping measure djj (cc, y) of the boundary process X. In particular, 
the Levy system characterizes the boundary process completely and we have obtained 
another probabilistic formulation of Calderon’s problem: Given pure jump processes on 
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dD generated by and Ak^, respectively. Show that the corresponding Levy systems 
coincide, i.e., -/V^i = if and only if = K2. 

5.2. Second version. An interesting aspect with regard to the first version is the fact 
that it may be translated into an equivalent parabolic unique continuation problem 
in terms of transition kernel densities of absorbing diffusion processes on D. More 
precisely, let ki,K 2 G be isotropic conductivities andp?; P2 the corresponding 

transition kernel densities. Define for any function cj) G the functions 

v{t,x)-.= f {p°{t,x,y)-p°2{t,x,ymy)dy, w{t,x) := j pl{t,x,y) 4 >{y) dy 
J D J D 

and set 

cit, x) := —V • {KiVwit, x)) + V • (/v2V?«(t, x)). 

By these definitions, we have 

{ dtv{t, x) = V • (KiV?;(t, x)) + cif, x), if, x) G (0, 00 ) x D 
u(0, x) = 0, x € D 

v{t,x) = 0, {t,x) G (0,oo) X dD. 

By the probabilistic formulation from above, the uniqueness result from Caro and 
Rogers [5] corresponds to the following assertion: Assume that 

poo 

(29) / d,.,,,(,,)d,,,,(y){p\[t,x,y)-pl[t,x,y))dt = d. 

Jo 

Then v must be identically zero for any (j) G L^{D). 

5.3. Third version. Let us conclude this section by deriving yet another equivalent 
formulation of Calderon’s problem which is based on the construction for the reflecting 
Brownian motion from [17]. To be precise, as in [17], we adopt Ito’s idea of regarding 
the excursions of X from the boundary as a point process taking values in the space 
of excursions. First, we recall some definitions from [17]. A measurable function 
e : [0, 00 ) —?► dD U {5} is called a point function if the set J(e) := {s > 0 : Cs G dD} is 
countable and we denote the set of point functions by 'P{dD). For each point function 
e, the counting measure on (0, 00 ) x dD is defined by 

(30) ne{B) := {{s,x) G B : Cs = x}, B C (0, 00 ) x dD. 

Given the probability space (D,JF^P,j,) with filtration {Ft > 0}, a function e : 
D —>■ V{dD) is called a point process if for each B C B{dD), the increasing process 
t !->■ ne((0,t] X B) is adapted to {Xt,t > 0}. 

Let denote the space of continuous excursions of X from a G dD to 6 G dD, 
i.e., the space of continuous paths e G C'([0, 00); D) such that e(0) = a and there exists 
an / > 0 such that e(t) G D for all 0 < t < 1 and e{f) = b for all t > 1 . If the space 
is equipped with the filtration a{e{s),e G LF“’^,0 < s <t), then the space of all 
excursions is given by 

(31) W= y LF“, 

where 

(32) IF“ := y 

b^OD ,b^a 
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The same proof as in [17, Proposition 4.4] shows that the random set of jump times 
7 ^ Xs} is a countable and dense set and that there is a constant c > 0, depending 
only on the domain H, such that after any given time t > 0, there are always infinitely 
many jumps of size at least c. Therefore, we set 

J := {s £ (0,oo) : t(s—) < t(s)} and l{s) := t(s) — t(s—) 


and define the point process of excursions of X as a hP-valued point process such that 


(33) 


es{t) 


^t+T{s — ) ; 
^t(s) ; 


t < l{s) 
t > l{s), 


if s G J, whereas Cg '■= d ii s ^ J. 

Clearly, the family {et,t < 0} is adapted to {Tt^t > 0}. The following definition 
yields a useful characterization of point processes. 


Definition 5.2. A cr-finite random measure he on the measurable space ((0,c») x 
dD,B{{0,oo)) X B{dD)) is called the compensating measure of a point process e if 
there is a sequence {BfcjfcgN C B{dD), exhausting dD, such that 

(i) the mapping t >->• he((0,t) x Bk) is continuous for every fc G N; 

(ii) ]Ehe((0,t) X Bk) < oo for every A: G N; 

(iii) for every B G B{dD) contained in some Bk, the process 


t !->• ne((0, t) X B) — he((0, t) X B) 
is a martingale with respect to {Bt,t >0}. 


The following theorem is a straightforward generalization of [15, Theorem 5.1] for 
the reflecting Brownian motion. It provides a compatibility condition, connecting the 
excursion law ¥a^b on which is completely determined by with the Levy 

system of X. 

Theorem 5.3. Let k : D ^ be a symmetric, uniformly bounded and uniformly 

elliptic conductivity with components Kij G i,j = 1, ...,d, such that k satisfies 

(Al). Then the point process of excursions of X admits the compensating measure 

(34) heiiO,t)xB)= [ {Bn{e : e(0) = X4}ds, 

Jo 

where the excursion law Qa from a G dD is a a-finite measure on W°‘ defined by 


(35) Q^{B} = f ¥a,b{Br\{e: e{l) = b}}N,,{a,b)da{b) 

JdD 

for every measurable B C 1T“. 

We have thus arrived at the following equivalent probabilistic formulation of Calderon’s 
problem: Given the boundary process X with its corresponding Levy system (Nk^{x, y) dcr(t/), id*), 
show that Pa,f, determined by p'^ is the unique excursion law such that (34) yields the 
compensating measure of the point process of excursions of some reflecting diffusion 
process on D. 
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6. Conclusion and outlook 

In this work, we have obtained a probabilistic formulation of Calderon’s inverse 
conductivity problem. This formulation comes in three equivalent versions and each 
of them may yield both a novel perspective as well as a novel set of (probabilistic) 
tools when it comes to studying questions related to the unique determinability of 
conductivities from boundary data. Indeed, it was shown in [17] that for the case 
K = 1/2 and under a certain consistency assumption, the reflecting Brownian motion 
can be reconstructed from its point process of excursions and the boundary process. 
However, the consistency assumption was derived by using the reflecting Brownian 
motion to begin with and as it is noted in [17], there might be other consistency 
assumptions leading to other constructions. Showing that there are no other consistent 
constructions is equivalent to the probabilistic inverse problem which will be the subject 
of future research. 


Acknowledgments 

The research of M. Simon was supported by the Deutsche Forschungsgemeinschaft 
(DFG) under grant HA 2I2I/8 -I 58306. The research of P. Piiroinen was supported 
by Academy of Finland (AF) under Finnish Centre of Excellence in Inverse Problems 
Research 2012-2017, decision number 250215. He has also been supported by an AF 
project, decision number 141075 and ERC grant 267700, InvProb. The authors are 
pleased to thank Martin Hanke for insightful comments on this work. 

References 

[1] M. Aizenman and B. Simon, Brownian Motion and Harnack inequality for Schrodinger Operators, 
Comm. Pure Appl. Math., 35 (1982), 209—273. 

[2] K. Astala, M. Lassas and L. Paivarinta, Calderon’s inverse problem for anisotropic conductivity 
in the plane, Comm. Partial Differential Equations, 30 (2005), 207-224. 

[3] K. Astala and L. Paivarinta, Calderon’s inverse conductivity problem in the plane., Ann. of Math. 
(2), 163 (2006), 265-299. 

[4] A. Bencherif-Madani and E. Pardoux, A probabilistic formula for a Poisson equation with Neu¬ 
mann boundary condition, Stoch. Anal. Appl., 27 (2009), 739—746. 

[5] P. Caro and K. Rogers, Global uniqueness for the Calderon problem with Lipschitz conductivities, 
Preprint, arXiv: 1411.8001 (2014). 

[6] Z. Q. Chen, M. Fukushima and J. Ying. Traces of symmetric Markov processes and their char¬ 
acterizations, Ann. Probab., 34, (2006), 1052-1102. 

[7] D. Dos Santos Ferreira, C. Kenig, M. Salo and G. Uhlmann, Limiting Carleman weights and 
anisotropic inverse problems. Invent. Math., 178 (2009), 119-171. 

[8] M. Fukushima. Dirichlet spaces and strong Markov processes, Trans. Amer. Math. Soc., 162 
(1971), 185-224. 

[9] M. Fukushima. On a decomposition of additive functionals in the strict sense for a symmetric 
Markov process, in Dirichlet forms and stochastic processes (Beijing, 1993), de Gruyter, Berlin 
(1995), 155-169. 

[10] M. Fukushima, Y. Oshima and M. Takeda. Dirichlet forms and symmetric Markov processes, 
Walter de Gruyter & Co., Berlin, 1994. 

[11] P. Grisvard, “Elliptic problems in nonsmooth domains”. Pitman, Boston, MA, 1985. 

[12] B. Haberman, Uniqueness in Calderons problem for conductivities with unbounded gradient. 
Preprint arXiv: 1410.2201 (2014). 

[13] B. Haberman and D. Tataru, Uniqueness in Calderon’s problem with Lipschitz conductivities, 
Duke Math. J., 162 (2013), 496-516. 

[14] M. Hanke, N. Hyvonen, S. Reusswig. Convex backscattering support in electric impedance tomog¬ 
raphy, Numer. Math., 117 (2011), 373—396. 


18 


Petteri Piiroinen and Martin Simon 


[15] P. Hsu. Probabilistic approach to the Neumann problem, Comm. Pure Appl. Math., 38 (1985), 
445-472. 

[16] P. Hsu. On the Poisson kernel for the Neumann problem of Schrodinger operators, J. London 
Math. Soc. (2), 36 (1987), 370-384. 

[17] P. Hsu, On excursions of reflecting Brownian motion, Trans. Amer. Math. Soc., 296 (1986), 
239-264. 

[18] I. Karatzas and S. E. Shreve. Brownian motion and stochastic calculus. Springer-Verlag, New 
York, 1991. 

[19] J.-P. Lepeltier and B. Marchal, Probleme des martingales et equations differentielles stochastiques 
associees a un operateur integro-differentiel, Ann. Inst. H. Poincare Sect. B., 12, (1976), 43-103. 

[20] J. Necas. Sur une methode pour resoudre les equations aux derivees partielles du type elliptique, 
voisine de la variationnelle, Ann. Scuola Norm. Sup. Pisa, 16 (1962), 305-326. 

[21] A. Pazy. Semigroups of linear operators and applications to partial differential equations. 
Springer-Verlag, New York, 1983. 

[22] P. Piiroinen and M. Simon. From Feynman-Kac Formulae to Numerical Stochastic Homogeniza¬ 
tion in Electrical Impedance Tomography Preprint, arXiv:1502.04353 (2014). 

Petteri Piiroinen, Department of Mathematics and Statistics, University of Helsinki, 

FI-00014 Helsinki, Finland 

E-mail address: petteri.piiroinen@helsinki.fi 

Martin Simon, Institute of Mathematics, Johannes Gutenberg University, 55099 Mainz, 

Germany 

E-mail address: simon@math.uni-inainz.de 


